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SIGNIFICANCE  AND  EXPLANATION 


Spontaneous  instability  of  laminar  flames  has  been  observed  in  many 
experiments.  A  primary  destabilizing  effect  is  known  to  be  thermal  expansion 
of  the  gas  passing  through  the  flame  front  and  causing  hydrodynamic 
instability  of  the  flame.  This  is  manifested  by  a  wrinkling  of  the  flame  and 
the  resulting  corrugated  flame  front  may  be  seen  to  propagate  at  a  constant 
velocity  with  its  wrinkled  shape  well  preserved.  The  speed  of  propagation  of 
the  wrinkled  flame  is  markedly  greater  than  that  of  the  undisturbed  flame. 

This  paper  demonstrates  the  above-described  phenomenon  through 
theoretical  means.  The  nonlinear  integro-differential  evolution  equation 
derived  by  Sivashinsky  for  spontaneous  hydrodynamic  instability  of  a  laminar 
flame  is  considered.  A  steady  progressive  long-wave  solution  of  that  equation 
is  shown  to  produce  a  wrinkled  flame  front  with  a  logarithmic  singularity  in 
slope  at  the  cusp  of  each  wrinkle.  The  calculated  increase  in  propagation 
speed  is  found  to  agree  well  with  experimental  observation. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sumswry  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A  CALCULATION  OF  WRINKLED  FLAMES 


H.  V.  McConnaughey,  G.  S.  S.  Ludford  and  G.  I.  Sivashinsky 
Introduction 

A  nonlinear  flame  stability  analysis  carried  out  by  Sivashinsky  (1977) 
considers  instability  of  a  laminar  plane  flame  to  disturbances  in  the  hydro¬ 
dynamic  field,  in  that  work,  a  nonlinear  integro-dif f erential  equation 
describing  the  evolution  of  the  perturbed  flame  front  is  derived.  The  present 
study  seeks  a  solution  of  that  equation  which  represents  the  experimentally 
observed  phenomenon  of  a  steadily  propagating  flame  with  corrugated  or 
"wrinkled"  profile  (Markstein,  1970;  Lind  and  Whitson,  1977;  Ivashchenko  and 
Rumyantsev,  1978;  Groff,  1982).  Attention  is  accordingly  focused  on  the 
simplified  equation 

F  _(n) 

2  ^on(n)  +  Y  — <*  n  -  V  ,  (1) 

h-h 

which  will  be  motivated  presently.  The  function  FQ(n)  is  the  shape  of  the 
disturbed  flame  front  which  has  been  assumed  to  be  a  steady  progressive  wave 
of  long  (but  arbitrary)  period;  V  is  the  constant  speed  of  propagation  of 
the  front;  n  is  the  space  coordinate.  These  quantities  are  all  dimension¬ 
less.  The  parameter  7  is  proportional  to 
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where  p  /p  is  the  ratio  of  the  density  of  the  burned  gas  to  that  of  the 
unburned  gas  (the  thermal  expansion  coefficient). 

It  is  found  that  equation  ( 1 )  does  indeed  describe  a  wrinkled  flame 
front.  The  purpose  of  this  work  is  to  show  that  the  slope  of  the  flame 
profile  has  a  logarithmic  singularity  at  the  cusp  of  each  wrinkle  and  to  give 
an  accurate  calculation  of  the  profile. 

Development  and  Discussion  of  Equations 

Equation  ( 1 )  is  obtained  from  the  more  general  equation 


P  _( n,  T) 
On 


fot  +  4Fonnnn  *  t  Fonn  +  2  Fon  +  Y  - 
s  n-n 


dn  -  0 


Here,  FQ(n,T)  is  the  perturbed  flame  front;  n  and  t  are  the  space  and 
time  variables;  is  a  positive  0(1)  constant  which  plays  an  unimportant 

role  in  the  present  discussion.  For  a  derivation  of  (2),  the  reader  is 
referred  to  Buckmaster  and  Ludford  (1982).  Equation  (2)  is  seen  to  be 
equivalent  to  the  aforementioned  equation  derived  by  Sivashinsky  (1977)  by 
recognizing  the  identity 


F  _( n,  T) 

|k|exp[ik(n-n)]F0(n,T)dndk  -  2  i£m  0 n ■_ —  an 

n-n 

(Buckmaster  and  Ludford,  1982). 

The  limit  Y  ♦  0  and  the  slow  scales  Yn  and  y2  t  are  now  considered. 

Consequently,  the  fourth  derivative  term  drops  out  of  (2).  Singling  out  the 

negative  coefficient  of  F„  _  then  yields 

onn 

f  _("n,  t) 

pot  “  7"  Fonn  +  2  Fon  +  Y  -  dn  "  0  •  (3) 

•  n-n 

Michelson  and  Sivashinsky  (1977)  approximate  the  solution  of  this  equation 
using  a  finite  difference  scheme .  Their  nuswrical  work  shows  that  an  imposed 
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initial  disturbance  of  large  characteristic  dimension  gives  rise  to  a  steady 
progressive  wrinkled  flame  front  of  the  same  characteristic  length. 

Information  about  the  steady  state  can  be  extracted  from  (3)  by 
analytical  considerations*  If  FQ(  n,  t)  is  taken  to  be  a  steady  progressive 
wave,  it  may  be  written 

F0<n,t)  -  -VT  +  F0(n>  , 

where  v  is  a  constant.  In  addition,  for  a  wave  of  long  period,  the  second 
derivative  term  may  be  neglected.  The  amission  of  PQ  ^  is  expected  to  be 
invalid,  however,  in  the  neighborhood  of  anticipated  cusps  in  the  flame 
profile  where  the  curvature  must  be  comparatively  large  in  order  to  effect  the 
sudden  change  in  slope. 

These  modifications  of  (3)  yield  equation  (1).  Next,  the  origin  is  fixed 
at  a  point  where  FQ|)  vanishes,  consequently  V  may  be  expressed  as 

f  _(n) 

V  -  Y  Cm  - -  <**  '  <4> 

h 

Substituting  (4)  into  (1)  gives  the  integral  equation 

J  y2(n)  ♦  £„  t*-  -  ;)y(h)d'n  -  o  ,  (5) 

n-n  n 

where 

y(n)  -  y"1  F0n(r>)  • 

A  periodic  solution  of  this  equation  is  sought  and,  without  loss  of 
generality,  the  period  is  taken  to  be  2w.  (a  solution  with  arbitrary 
period  2p  can  then  be  obtained  by  changing  the  argument  from  n  to  pV*i 
the  integral  in  (4)  shows  that  V  is  unaltered.  Hence,  all  periodic  profiles 
propagats  at  the  same  speed  and  have  the  same  amplitude}  only  their  wavelenths 


are  different.) 


Equation  (S)  may  be  rewritten  as 


\  y2"»  -  -  I  tJ— - 


when  y(n)  has  period  2*,  this  is  equivalent  to 


\  y2(n)  -  -  I  +  T?ki)yU)dC 

n-~® 

“  "  2  Icot(^T3)  “  cot  -jlyC C)dC  .  (6) 

The  principal-value  integral  in  (6)  has  singularities  at  n  ■  ±w,  where 

y2(n)  ~  cot(^jp)y(  £)d£  •  (7) 

The  singular  behavior  there  is  found  to  be 

y(n)  ~  2  sgn(n?v)  ln|n?*|  near  n  “  *»  .  (8) 

To  see  this  at  n  -  *#  first  note  that 

Hi  cot^lytOd?  -  tf2*  cot(-^)y(C)d5 

~  y{  ^)d*  near  n  “  5  “  *  • 

It  is  easily  found  that 

^■.^■.^,.,,.1^  n«ar  «-  .  . 

hence  relation  (8)  satisfies  relation  (7)  for  n  ■  *»  the  result  for  n*  -» 
follows  similarly.  The  function  y( n)  is  thus  seen  to  have  a  logarithmic 
singularity  at  odd  factors  of  i,  i. e. ,  at  the  cusps  of  the  flame  front. 

The  singular  nature  at  the  cusps  of  the  flame  profile  is  reflected  by  the 
pertinent  numerical  results  of  Michelson  and  Sivashinsky  (1977).  They  took 
the  term  FQj1fl  (see  equation  (3))  into  consideration  in  their  calculations 
and  found  the  curvature  of  the  flame  to  be  surprisingly  high  at  the 
wrinkles.  The  singularity  there  creates  a  stagnation  zone  near  each  cusp  in 


the  burned  gas  region,  as  occurs  at  the  tip  of  a  Bunsen  flame  (Figures 
1a,b).  The  structure  of  the  flame  front  at  a  wrinkle  (Figure  la),  however,  is 
quite  different  from  that  of  a  Bunsen  wedge  flame  (Figure  1b). 

Assume  now  that  the  wrinkled  flame  profile  is  symmetric,  i.e.,  that 
y(D)  is  an  odd  function.  Equation  (6)  may  then  be  written 

y2(n)  m  ~  (cot(4“b  +  cot(-^j^)  -  2  cot  ^]y(£)d£  •  (9) 

This  equation  describes  an  odd  and  2*-periodlc  solution  of  equation  (1).  It 
is  a  simple  exercise  to  show  that  equation  (9)  therefore  represents  the 
solution  of  the  problem  addressed  by  Sivashinsky  (1976)  in  his  discussion  of 
corrugated  flame  fronts,  in  that  analysis,  the  following  equation  for  the 
shape  of  a  two-dimensional,  long-wave  corrugated  profile  is  derived 

m 

A  -  G2(n)  ■  “jjj  J*w  G(C)sin  nC  cos  nndC  .  (10) 

n-1 

Here,  the  function  G(n),  which  is  proportional  to  the  slope  of  the  flame 

front,  is  also  assumed  to  be  2»-periodic  and  odd.  The  constant  A  measures 

the  speed  of  propagation  of  the  front.  This  equation  is  recovered  from  (1)  by 

writing  F  _( n)  as  the  Fourier  sine  series 
on 

m 

“  I  sin  nn  F  (Osin  n?dC  , 

i*«1  v 

so  that  the  integral  in  (1)  nay  be  written 

l  llw  F  (Oain  nCdC  £  ^  .  (11) 

n-i  n-n 

The  principal-value  integral  in  (11)  is  equal  to  v  cos  no.  Hence,  equation 
(10)  is  obtained  by  defining  G(n)  ■  F0r)(n)/2irf  and  A«  V/(2rf)2. 

In  the  analysis  of  Sivashinsky  (1976),  the  solution  of  (10)  is  crudely 
approximated  by  the  solution  of 

A  -  j  G2(H)  ■  /*w  6(  O  sin  C  cos  ndc  . 
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In  the  present  work#  an  accurate  solution  is  found  by  a  numerical  calculation 
of  the  solution  of  (9)  in  which  the  appropriate  singular  behavior  at  v  is 
imposed. 

Numerical  Solution 

From  equation  (9)  and  result  (8)#  a  system  of  N  ♦  1  equations  is 
obtained  for  y^  “  y(x^),  where  x^  ”  (i-Dx^/N,  i  ■  1,2,...,N+1  and  the  end 
value  xf  is  less  than  but  close  to  v.  The  principal-value  integrals  in  (9) 
are  approximated  near  their  singularities  by  expanding  y(n)  in  a  Taylor 
series,  which  produces  ordinary  integrals.  All  integrals  are  then 
approximated  by  simple  quadrature  and  all  derivatives  by  finite  differences. 

In  this  way,  the  following  system  (12)  is  obtained.  This  system  of  nonlinear 
algebraic  equations  is  solved  numerically  using  the  Newton- Raphson  method. 


x.  i-2  N  ( j-i)x.  N  (j+i-2)x 

«  l  *  1  Hy^ootl— 5,  n  +  E  y.i(cotl — 2m — ~ 

"  j-2  J-i+2  3  j-2  3 

(j-Dx.  y_. .  (*-i+1)x  (M+i-Dx 

2  cott  -a|  —  1)  ♦  (cot l - — - -]  ♦  cot[ - — - 


2  cot  /)}  ♦  25’  -  **2  •  ° 


(12) 


(i  ■ 

*«  ■  -»  u<-v 

Vi  ■  -2  • 

The  last  two  equations  in  (12)  represent  the  logarithmic  singularity  at  *. 

For  this  calculation,  N  was  taken  to  be  27  and  x£  -  3.13.  imposing 
the  behavior  (8)  at  x^  alone  was  found  to  yield  unacceptable  results?  this 
occurred  for  all  three  values  of  M  that  were  tested  (N  -  9,  18,  27). 

~7- 


Fixing  y  at  the  last  two  points  according  to  (8)  corrected  this,  however 
Fixing  the  third-to-last  value  as  well  did  not  appreciably  change  the 
solution,  which  is  shown  in  Figure  2.  Interpolating  between  the  points 


*xi'yi*  an<J  numerically  integrating  from  0  to  x£  gives  the  approximate 
shape  of  the  wave  front 

-  FQ(n)  -  /q  y(x)dx  ,  ne  [o,t)  , 

as  illustrated  in  Figure  3.  A  smooth  curve  is  fitted  to  the  points 
(x^,  y  1  Fq ( xi ) >  end  is  then  appropriately  extended  beyond  the  interval 
[0,x£]  to  generate  the  graph  in  Figure  4  of  a  portion  of  the  complete  flame 
profile. 

Discussion  of  Results 

The  constant  speed  of  propagation 

V  -  Y 2  J*  cot  ^  y(x)dx 

may  be  computed  from  the  numerical  results  obtained  for  y(x).  It  is  found 
w  2 

that  V  -  1.4Y  .  The  corresponding  dimensional  speed  of  the  wrir'xled  flame 
may  be  shown  to  be 


0 

w 


°b{1  ♦  ')2—2 

(2»r 


where  a  is  the  speed  of  the  undisturbed  flame  front  relative  to  the  burned 

gas.  A  typical  value  of  the  thermal  expansion  coefficient  p,_/p  is  0.2. 

b  u 

Thus,  the  speed  of  the  steady  plane  flame  is  amplified  by  a  multiplicative 
factor  of  1.6.  This  is  in  excellent  agreement  with  recent  experimental 
observations  of  large-scale  flames  propagating  in  thin  spherical  or 
hemispherical  shells  (Ivashchenko  and  Rumyantsev,  1978;  Lind  and  Nhitson, 
1977).  In  these  investigations,  it  was  observed  that  for  mixtures  with 
markedly  different  burning  velocities,  the  measured  space  velocity  was  1.5 
to  2  times  the  normal  burning  velocity  measured  in  the  laboratory.  Also, 
the  flame  became  rough  as  it  expanded,  having  a  "pebbled”  appearance. 


Numerical  investigation  by  Michelson  and  Sivashinsky  (1977)  of  equation 
(3)/  which  includes  curvature  effects,  yielded  the  result 

U  -  U.  {1  +  (1  -  ~)2(.18) }  . 

w  b  p 

u 

For  p^/ »  .2,  this  produces  a  magnification  factor  of  only  1.1.  This  is 

apparently  due  to  the  narrowness  of  the  integration  interval  taken.  Over  a 

relatively  short  interval,  the  stabilizing  effect  of  the  curvature  term  F_ 

0  nh 

may  be  significant,  and  so  the  wrinkles  may  be  underdeveloped.  Indeed,  it  is 
known  from  experiments  on  spherically  expanding  flames  that  the  speed  of  a 
wrinkled  flame  stabilizes  at  an  extremely  low  rate  (Palm-Leis  and  strehlow, 
1969). 

Another  system  to  which  the  present  findings  may  have  some  relation  is 
flame  propagation  between  two  parallel  walls  (Uberoi,  1959).  In  particular, 
such  a  flame  may  be  described  by  the  curved  front  found  here  between  two 
consecutive  cusps.  *~he  propagation  of  a  flame  in  a  channel  was  recently 
considered  by  zeldovich  et  al.  (1980)  where,  however,  estimates  obtained  both 
for  the  flame  shape  and  for  its  propagation  velocity  were  comparatively 
coarse.  In  that  work,  stagnation  zones  near  the  channel  walls  were  taken  into 
account,  as  shorn  in  Figure  5. 
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